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Hard X-ray emission from a young massive star-forming cluster 
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Abstract. We report the detection of hard X-ray emission (> 2 keV) from a number of point sources associated 
with the very young massive star-forming region IRAS 194104-2336. The X-ray emission is detected from several 
sources located around the central and most deeply embedded mm continuum source, which remains undetected 
in the X-ray regime. All X-ray sources have K-band counterparts, and those likely belonging to the evolving 
meissive cluster show near-infrared colors in the 2MASS data indicative of pre-main-sequence stages. The X-ray 
luminosities around 10^^ ergs~^ are at the upper end of luminosities known for low-mass pre-main-sequence 
sources, and mass estimates based on the infrared data indicate that at least some of the X-ray detected sources 
are intermediate-mass objects. Therefore, we conclude that the X-ray emission is due to intermediate-mass pre- 
main-sequence Herbig Ae/Be stars or their precursors. The emission process is possibly due to magnetic star-disk 
interaction as proposed for their low-mass counterparts. 
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1. Introduction 

In the past, star formation research in the X-ray regime 
has fo' 

stars, 'pliuse ubjetLs uuuL mainly in Llie sufL 



used sLiungly uu luw-mass ubjetLs, e.g., T Tauii 



y in Llie sulL lange uf the 
with typical X-iay lumino&i- 



Furthermore, X-ray variability is observed in all types of 
low-mass pre-main-sequence objects. For an excellent re- 
cent r eview on t hese topics see Feigelson &: Montmerk 
(1999D. Recently, [Tsuboi et al. (200l|) and [Tsujimoto et 



X-ray spectmm (■- 
ties between 10^* — 10"^" ergs^^. The satellite observatory 
ROSAT was an ideal instrument to study such stars within 
a few 100 pc distance from the Sun, and the observed 
emission can be explained by enhanced sola r-type mag- 
netic activity ( Feigelson fc Montmerlc, 1999 ). In the last 
few years, a rising number of Class I protostars, which 
are still deeply embedded within their natal molecular 
cores {Ay « 10 — 100 mag), have been detected in the 
hard X-ray regime between 2 and 10 keV with X-ray lu- 
minosities higher than 10^° ergs~^. These detections were 
mostly made with the X-ray satellites ASCA and, most re- 
cently, with Chandra. Magnetic star-disk interactions are 
thought to b e the most likely explanation for the hard X - 
ray emission (iHayashi et al., 1996| ; pvIontmerle et al., 2000| ). 



al. (2002) reported the first tentative detections of deeply 



Send offprint requests to: H. Beuther 

* Based on observations with the Chandra X-ray observa- 
tory, the Calar Alto 3.5 m telescope and the IRAM Plateau de 
Bure Interferometer. Chandra is operated by the NASA and the 
CXC, Calar Alto by the Max-Planck-Institut for Astronomy in 
Heidelberg, jointly with the Spanish National Commission for 
Astronomy, and IRAM is supported by INSU/CNRS (France, 
MPG (Germany), and IGN (Spain). 



embedded Class protostellar candidates in 0MC3 by 
Chandra. 

In comparison with the low-mass regime, X-ray obser- 
vations of massive star-forming regions have been rare. 
Due to the high visual extinction within such regions 
{Ay up to a few 100 or even 1000), soft X-ray emis- 
sion is completely absorbed by the gas along the line of 
sight. iHofncr fc ChurchweU (1997[ ) detected with ASCA 
for the first time hard X-ray emission in the massive star- 
forming region W3. Because of the low angular resolution 
of ASCA (> 1') they could not determine whether the 
emission is caused by the superposition of many point 
sources, e.g., protostellar clusters, or whether it is due 
to a wind-shocked cavity resulting from strong stellar 
winds interacting with the surrounding medium. Recently, 
[ChurchweU (2001 ) reported that Chandra data of the same 
region with a spatial resolution of 0.5" resolve the emis- 
sion into many individual sources distributed over the en- 



tire W3 complex. Similarly, Garmire et al. (2000) and 



Feigelson et al. (2002) reported about 1000 X-ray emit- 
ting pre-main-sequence stars between 0.05 and 50 Mq 



in the Orion Nebula. Additionally, Zinnecker & Preibisch 
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(1994) found X-ray emission with the ROSAT satellite in 
the soft X-ray band associated with several intermediate- 
mass Herbig Ae/Be pre- main-sequence stars. The derived 
X-ray luminosities for the Herbig Ae/Be stars are rang- 



high for the far kinematic distance (Beuther et al., 2002b). 
At the distance of 2.1 kpc, its infrared derived luminos- 
ity is 10'* L0, and one observes two adjacent star- forming 
cores with masses of 840 Mp, and 190 Mf^ ( ^ridhara n et 



ing be;wecn 10'^" ergs * and lO'^'' ergs ^. Preibisch & al., 2002; Beuther et al., 2002a). Each core drives a mas 



Zinneccer (1995) speculated that the emission might orig- sive bipolar outflow in east- west direction (Beuther et al 



inate f'om coronal activity due to shear dynamo action. 



A study of the more distant molecular clouds Monoceros 
and Rosette also found indirect evidence for X-ray emis- 
sion from intermediate-mass pre-main-sequence sources 
( Grcgorio-Hctcni ct al., 199^ ) . An example of hard X-ray 
emission from a Herbig Be star is MWC 297 (Hamaguchi 
et al., 2000). In a recent X-ray study of the Monoceros 



R2 molecular cloud, Kohno et al. (2002 ) detected X-ray 
emission for stars of all masses, in particular hard X-ray 
emission from high-mass pre-main-sequence or Zero-age- 
main-sequence stars. Based on these results. X-ray emis- 
sion seems likely to be an ubiquitous phenomenon in the 
protostellar evolution of stars of all masses. 

So far. X-ray studies in high-mass star-forming re- 
gions focused on more evolved massive star formation 
sites, while the very deeply embedded phase - and thus 
the youngest stage of stellar evolution - has not been 
detected at all. Results obtained in recent years by our 
group ( ISridharan et al., 2002| ; [Beuther ct al., 2002a| ,[bp|) 
and other groups studying massive star formation (e.g 
Cesaroni et al. 1997| , [Zhang et al. 2002| , [Tan fc McKee 
2002 |, ^oike 2002| ) support the hypothesis that massive 
stars form via disk accretion in a similar fashion as low- 
mass stars. Therefore, high-mass star-forming cores are 
promising candidate regions where in rather small spatial 
areas a number of sources could be hard X-ray emitters 



via the physical process of star-disk interactions (Hayashi 
et al., 1996; Montmerle et al., 2000| ). Necessary obser 



vational requirements are first of all sensitivity in the 
hard X-ray regime (> 2 keV), because only hard X- 
ray photons can penetrate high gas column densities. 
Additionally, high angular resolution is needed to resolve 
different sub-sources of the forming cluster. While ROSAT 
was not sensitive to hard X-ray photons, the spatial reso- 
lution of ASCA was not sufficient to study massive star- 
forming regions in detail at their typical distances of a 
few kpc. The new-generation X-ray satellite telescopes 
Chandra and XMM-Newton comprise both features, be- 
ing sensitive up to 10 keV, and having a spatial resolution 
of 0.5" (Chandra) and 15" (XMM-Newton), respectively. 
Especially Chandra is able to resolve many different sub- 
sources as impressively demonstrated in Orion by Garmire 
et al. ( ^000[ ) and [Feigelson et al. (2002| ). 

Here, we present a Chandra X-ray study of the very 
young, massive and deeply embedded star-forming clus- 
ter IRAS 19410-f2336. IRAS 19410-1-2336 is part of a 
large sample of 69 high-mass protostellar candidates which 
has been studied extensively in a series of papers dur- 



ing the last years (Sridharan et al., 2002; Beuther et al 
2002a[]p[[c[,p[). We assume the source to be located at its 



2002b). At the center of the southern massive core, a very 
compact and weak ('^ 1 mJy) cm wavelength source is 
detected, which coincides with H2O and Class 11 CH3OH 



near kinematic distance of 2.1 kpc ( Sridharan et al., 2002| ), 
because the derived outflow parameters are unreasonably 



maser emission ( [Beuther et al., 2002c[ ). Figure 0(top) gives 
an overview of the region of interest with the 1.2 mm 
dust continuum data ( Beuther ct al., 2002a[ ) superposed 
on an infrared K-band image (§2.3). As the source is lo- 
cated in the Galactic plane, confusion due to foreground 
and background sources is expected and has to be dis- 
entangled by the different observations. We focus on the 
X-ray emission of this region and correlate the detected 
X-ray sources with high-resolution images in the mm and 
near-infrared regime. §2 describes the different observa- 
tions we performed (X-ray, near-infrared and mm data), 
and in §3 we derive the main physical parameters of this 
star-forming region from our data (source detections, spec- 
tra. X-ray luminosities, plasma temperatures and masses). 
Finally, §4 compiles our conclusions and puts the observa- 
tional findings into a more general framework of high-mass 
star formation. 

2. Observations 

2.1. Chandra X-ray observations and data reduction 

IRAS 19410-^2336 was observed with the ACIS-S3 chip on 
board the X-ray telescope Chandra for 20 ksec on October 
15 2001. The ACIS-S3 aim-point was centered at the IRAS 
position R.A. 19:43:11.4, Dec. 23:44:06.0. (J2000.0) that 
is coincident with the main mm emission peak. The field 
of interest shown in Fig. |l| is much smaller than the ^ 
8' X 8' field of view of the S3 chip, and the FWHM of 
the point-spread-function of this region is approximately 
0.5". The data were reduced with the CIAO 2.2 software 
package using the CALDB 2.10 database. Both packages 
are provided by the Chandra X-ray center]^. The basic data 
product of our observation is the level 2 processed event 
list provided by the pipeline processing at the Chandra 
X-ray center. 

Light curves were extracted for the whole S3 chip as 
well as just for the region presented in Figure |l|. No flares 
or enhancements due to low energy protons are observed. 
The whole observation time is usable as good time in- 
terval. The background emission is about 0.1 cts pixel^^, 
but as one or two photons could be detected by back- 
ground fluctuations anyway, we conservatively use 6 cts as 
the detection threshold for an X-ray source in this field. 
A wavelet based source detection algorithm was applied 
using default parameters (Freeman et al., 2002), and 13 
sources between 7 and 19 cts are found in the field of in- 
terest (Table |^). Based on the 1.2 mm continuum map 

^ Mtp:/ /asc. harvard.edu/ciao/download. html 
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sequence sources (|Feigelson fc Montmerle, 1999| ; Preibisch 



R.A. [J2000.0] 



12 13 



■■ ■■■'*„■, 



11 

*io > 



*9 



©2 




19 43 15 10 

R.A. [J2000.0] 

Fig. 1. top: The contours present the large-scale 1.2 mm dust 
continuum emission from 10% to 90% (steps of 10%) of the 
peak flux ( Beuther et al., 2002e ). The grey-scale shows the 
K-band image, bottom: The contours show the PdBI 2.6 mm 
continuum dust cores (black: levels 3.6(1.2=lcr)8.4 mjy/beam, 
white: 9.6(2.4)20 mJy/beam) superposed on the K-band im- 
age. Circles and stars mark the X-ray sources presented in this 
paper. Sources with circles are most likely associated with the 
massive star-forming region, whereas the asterisks represent 
those that might be foreground objects. The box outlines the 
region presented in Fig. W. 



from the 30 m telescope ~ assuming optically thin dust 
emission -, we can estimate an average hydrogen column 
density of iVn = 5 x 10^^ cm~^ of the molecular cloud (for 
details on the dust emission see Beuther et al. 2002a). A 



Raymond-Smith plasma model (Raymond & Smith, 1977) 
with a 3 keV plasma, which is commonly used to model 
the X-ray emission of very young and embedded pre-main- 



& Zinncckcr, 2002), 6 cts correspond to an unabsorbed flux 
limit of 2 X 10^^* ergcm^^s^^. Based on the logA^-logS 
distribution of Giacconi et al. (2001), in a field as small 



as ours less than 0.1 extragalactic background sources are 
expected. The source positions in the field of interest are 
compared with 2MASS near-infrared data and should be 
correct within less than 1". 



2.2. Plateau de Bure mm observations 

We observed IRAS 19410-1-2336 in Summer 2001 with the 
Plateau de Bure Interferometer (PdBI) at 2.6 mm in the 
D (w ith 4 antennas) and C (with 5 antennas) configura- 
tion ( Guilloteau ct al., 1992 ). The simultaneously observed 
1 mm data were only used for phase corrections because 
of the poor Summer weather conditions. The 3 mm re- 
ceivers were tuned to 115.27 GHz (USB) (centered at the 
^^CO 1^0 line) with a sideband rejection of about 5 
dB. At this frequency, the typical SSB system temper- 
ature is 300 to 400 K, and the phase noise was below 
30°. Atmospheric phase correction based on the 1.3 mm 
total power was applied. For continuum measurements, 
we placed two 320 MHz correlator units in the band to 
cover the largest possible bandwidths. In this paper we 
are focusing on the 2.6 mm continuum data, the CO 
line observations will be presented elsewhere. Temporal 
fiuctuations of the amplitude and phase were calibrated 
with frequent observations of the quasars 1923-1-210 and 
20234-336. The amplitude scale was derived from measure- 
ments of MWC349 and CRL618, and we estimate the final 
fiux density accuracy to be ~ 15%. To cover both cores 
a mosaic of 10 fields was observed. The final beam size is 
4.1" X 3.6" (P. A. -108°). We obtain a 3cr rms of - 5 mJy, 
which corresponds to a mass sensitivity limit of approx- 
imately 15 Mq as suming optically thin dust emission at 



2.6 mm (see, e.g., Beuther ct al. 2002a for the deviation 



of dust parameters from mm observations). 



2.3. Calar Alto near-infrared observations 

The near-infrared camera Omega Prime on the 3.5 m tele- 
scope on Calar Alto/Spain was used to obtain K' wide field 
images of the region around IRAS 19410+2336 in June 
2001. At a pixel scale of 0.4" pixer^, the 1024x1024 pixel 
array provides a field-of-view of ^ 6.8' x 6.8'. A five posi- 
tion dither pattern was applied to image the field and to 
allow for a correction of array defects that were identified 
from well illuminated flat-field frames (for the dark pix- 
els) and a dark frame (for the hot pixels). In each dither 
position, 15 individual exposures of 2 sees were stacked 
into one frame of 30 sees total integration time. For each 
position, the thermal background (sky) was computed by 
median combining the four frames resulting from the other 
dither positions. This sky frame was subtracted from the 
respective science frame, which also removes the bias level. 
The resulting frame was divided by a flat-field resulting 
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from dome flats (the difference of a number of frames of 
the dome illuminated by a tungsten lamp and a number 
of frames taken without illumination). The images of the 
various dither positions were then registered and averaged 
into the final image. The total integration time in the cen- 
ter of the final image is 2.5 min, which yields a detection 
limit (5(7 peak flux) for point sources of roughly K'=17.5. 
The seeing-limited angular resolution in the final image 
is about 1.2". Because the observation was conducted at 
high air-masses, its photometry is rather poor. Thus, we 
used the Omega Prime data for identification only, and 
the photometry is taken from the 2MASS catalog. 



3. Results 

3.1. Source detections and identification 

Fig. |l|(bottom) shows that in the very young massive star- 
forming region IRAS 19410-1-2336 hard X-ray emission 
from several point sources is detected. No evidence for 
extended or diffuse emission is found within our detection 
hmits (§|]|). At mm wavelengths, the two massive cores 
detected with the IRAM 30 m single-dish observatio ns at 
an angular resolution of 11" (Beuther et al., 2002a) split 
up into many sub-sources when observed at higher spa- 
tial resolution (Fig. ^ bottom), confirming that a clus- 
ter of stars is forming. The interferometric data highlight 
the real massive cores, whereas the large-scale surround- 
ing core emission is filtered out by the interferometric ob- 
serving technique. A detailed analysis of the protocluster 
properties will be presented elsewhere. Here we are focus- 
ing on the spatial associations of X-ray, near-infrared and 
mm sources. 

Thirteen X-ray sources are detected within the field of 
interest (Table |l|). Remarkably, nearly all X-ray sources 
are, within the positional uncertainty, associated with 
near-infrared counterparts observed in the K-band (X3 
and X5 only tentatively). Because the K-band extinction 
is approximately only 10% of the visual extinction, in 



the characteristics of the associated near-infrared coun- 
terparts (Fig. H). While a few of the near- infrared sources 
are located around the unreddened main sequence in the 
color-color diagram, many sources are reddened and show 
near-infrared excess, suggesting that they are pre-main- 
sequence objects surrounded by circumstellar material. 
This is expected in view of the early evolutionary stage, 
the large core masses and the high column densities of the 
region (Beuther et al., 2002a). Table |l| compiles reddening 
estimates of the X-ray near-infrared counterparts. The 
of XI is most uncertain because XI has the largest near- 
infrared excess which makes the determination more 
difficult. Using the reddening scale (Fig. we can still 
estimate an approximate value of the visual extinction Ay 
of XI. 




H-K 



Fig. 2. Near-infrared color-color diagram of X-ray associated 
sources based on 2MASS data. The numbers correspond to the 
sources numbered also in Fig. hi. The full line present s main 



the K-|band we can observe more embedded regions of sequence colors with three labeled spectral star types ( pucat: 



star-forming cluster than possible in the optical, We rc- 



call that the extinction in the soft X-ray regime is sim- 
ilar to the K-band extinction, and that it drops even 



et al., 2001), and the d ashed lines show the redde ning band for 
main sequence colors (Rieke & Lebofsky, 1985). The dotted 



further going to higher X-ray energies (Casanova et al 



1995 ; Ryter, 1996 ). In contrast to the K-band counter- 



parts, only 2 sources (X2 and X8) are found in the direct 
vicinity of mm dust cores that are tracing the coldest and 
deepest embedded regions where presumely the youngest 
sources are found. The whole region is also covered by the 
2 Micron AU Sky Survey (2MASS0) in the near-infrared 
bands J, H and K. Except for sources X2, X3, X5 and 
XI 1, which are detected by our K-band Omega Prime ob- 
servation but which are too faint to be detected by the 
2MASS survey, all other X-ray sources are also 2MASS 
point sources. A color-color diagram helps to determine 



For details 



the 2MASS web-site 



at 



http:/ /www. ipac.caltech.edu/2mass/releases/docs. html 



line gives the reddening scale with encircled visual extinctions 
([mag]). Source XI shows the strongest near-infrared excess. 



3.2. Spectral information of the X-ray data 

3.2.1. X-ray spectra and X-ray luminosities 

The number of photons per detected source ranges from 7 
to 19 counts within the 20 ksec observation time (Table |l]). 
We note that we are dealing with low- number statistics, 
and this has to be taken into account when deriving the 
physical parameters of this region. Figure ^ presents the 
photon energies of the 13 X-ray sources versus time. As 
expected, we detect emission mainly between 2 keV and 
6 keV in the near-infrared reddened and thus embedded 
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sources, whereas the sources, that are near main sequence 
colors in the color-color diagram, emit more in the soft 
regime below 2 keV (see also mean values in Figure |3| and 
Table p. The deficiency of soft photons (< 2 keV) in the 
direction of the embedded sources is most probably due to 
the high column densities that absorb the soft X-rays. The 
lack of photons above 6 keV may be a real feature of the X- 
ray sources, but it has also to be noted that the sensitivity 
of Chandra above 6 keV decreases significantly. The effec- 
tive collection area of Chandra at 8 keV drops to less than 
a quarter of that at 4 keV|^ The fact that we detect hard 
X-ray emission (> 3 keV) of the embedded sources fits in 
the general finding that the younger the sources are, the 



harder the X-ray emission is (e.g., [Feigelson fc Montmcrle 
199g[o} §0). The X-ray spectra of XI to X7 show all such 
hard X-ray signatures, which makes them likely to belong 
to the forming cluster. The distribution of photon ener- 
gies cannot be explained by a low-temperature- Raymond- 
Smith plasma of 1-2 keV ( Raymond fc Smith, 1977 ), usu- 



ally used to fit the X-ray spectr a of low-mass T Tauri stars 
( Feigelson fc Montmerle, 1999 ), because the column den- 
sities A^H, that are necessary to attenuate all the emission 
below 2 keV(A^H a few times 10^'^ cm^^), would imply av- 
erage X-ray luminosities as high as 10^'' ergs~^. This value 
exceeds the known luminosities of young stellar objects by 
orders of magnitude (§0), and it could hardly be reached 
even during strongest flaring events (e.g., IPreibisch et al 



1993 |, |(prosso et al. 1997| , [Tsuboi et al. 20001 ). Contrasting 
these hard spectra, X8 to X13 show softer spectra. While 
some of them might be unrelated objects, others could 
belong to the cluster but be in a more advanced state of 
evolution (for details see §3.4). 

In order to derive estimates of the X-ray plasma tem- 
peratures and the X-ray luminosities, we performed a de- 
tailed spectral analysis. We extracted individual pulse- 
height spectra and built the corresponding redistribution 
matrices and ancillary response files for each source. We 
fitted the ungrouped pulse- height spectra with XSPEC, 
using a single-temperature thermal plasma model (model 
"raymond", Raymond & Smith 1977) plus the absorption 



model "wabs" . The plasma temperatures and column den- 
sities derived in these fits are compiled in Table |^. The 
fitting results were then used to integrate (in XSPEC) the 
model fiuxes over the 0.2-10 keV band and to compute 
the X-ray luminosities. In some cases the low number of 
detected photons allowed only to derive lower limits to 
the plasma temperature. We are fully aware that, given 
the rather low number of photons per source, the derived 
spectral parameters inherit a relatively large statistical 
uncertainty. Nevertheless, we believe that the fitting re- 



(kT > 3keV) and high colum densities (A^h > 10^^ cm^^); 
in contrast, the sources X8-X13 show low plasma temper- 
atures (kr < 1.5 keV) and (with the exception of X8) 
low colum densities {Nu < lO^^cm"^). These two groups 
exactly reflect the spatial distribution of the sources: XI- 
X7 are located inside the central region of the molecular 
clump, while X8-X13 lie outside the molecular clump (Fig. 

0). 

(2) The column densities derived in the spectral fits 
agree in nearly all cases rather well with expectations from 
the infrared colors of the sources (Table I, Av = iVn/Z x 
10^1). 

(3) iGetman et al. (2002| ) have performed fitting simu- 
lations of pulse-height spectra with low numbers of counts 
to estimate the statistical uncertainties of the fitting re- 
sults. They found that even if the column densities and 
temperatures derived in the fits are subject to rather large 
uncertainties, the broad-band X-ray luminosities derived 
from the fitting results are quite reliable (statistical error 
< 35% for spectra with 15 counts). 

Very likely, the X-ray emission of these objects origi- 
nates from a hot thermal plasma typical for very young 
embedded objects (e.g., p^eigelson fc Montmerle 1999 



Preibisch fc Zinnecker 2002| , |Kohno et al. 2002| ). But as 



the number of photons and by that the spectroscopic in- 
formation is small, other X-ray emission processes have to 
be taken into account, e.g., non-thermal power-law distri- 
butions as a signature of synchrotron emission, which has 
been observed in the radio regime by Reid et al. (1995 ) 
in the massive star-forming region W3(0II). Therefore, 
we also estimate X-ray luminosities with three different 
power-law distributions oc E~°' {a — 0.5, 1, 2) to derive 
an estimate for the influence of a on the total X-ray lumi- 
nosity (using the web-based PIMMS toolQ). The derived 
luminosities vary by less than a factor 1.5 for the differ- 
ent a. As hot plasmas are dominated by bremsstrahlung 
above kT - 1 keV ( |Fink fc Triimper, 198^ ), which can 
also be approximated by non-thermal power-law distribu- 
tions, the X-ray luminosities derived this way are close 
to the luminosities calculated using the Raymond-Smith 
plasma (they differ only a factor 2 at maximum). 

We calculated only for XI to X8 (Table H), because 
the other sources could be foreground or background ob- 
jects, and thus no reliable distance is known. We find X- 
ray luminosities in the 10'^^ ergs~^ range. These values 
are high, and above or at the upper end of classical T 
Tauri X-ray luminosities, i.e., in the regime of younger 
Class I sources ( Feigelson fc Montmerle, 19991) and/or 



suits cdJiiLaiii valuable liifuimaLluu, fui Ihu ruUuwiiig Ihn 



reason! r 



more massive Herbig Ac/Be type objects (Zinnecker fc 



(1) The spectral parameters derived ni rhc firs al- 
low a clear distinction of two different groups among the 
sources: sources X1-X7 show high plasma temperatures 



Preibisch, 1994; Gregorio-Hetem et al., 1998 



Hamaguch: 



et al., 200C). Unresolved binary systems could alter the 



results slightly, but not by orders of magnitude. 



^ See Chandra Pro posers' Observatory Guide, 2000, 

http:/ /cxc. harvard.edu ^ http:/ /asc. harvard.edu/tooIkit/pimms.jsp 
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Fig. 3. Photon energy versus time for the 13 detected sources within the field of interest. The dotted hues mark the mean 
detected photon energy in each source (see also Table m) . 



Table 1. Compilation of the parameters derived from the X-ray and 2MASS data: positions of the X-ray sources, number of 
counts per source, mean energy of detected photons, best fits of Raymond-Smith plasma temperatures (kT [keV] and T [K]) and 
derived X-ray luminosities Lx with error ranges. Colons mark less conclusive fits. Additionally, the column densities A^h (and 
error ranges) from the X-ray fits and the 2MASS derived extinctions are presented. 



source 


R.A. 


Dec. 


# 


mean 


kT 


T 


Lx 


A^H 






(J2000.0) 


(J2000.0) 


[cts] 


[keV] 


[keV] 


[10^ K] 


[lO^^ergs-^] 


[lO^^cm-^] 




XI 


19:43:10.8 


-f 23:44:04. 6 


17 


2.6 


5.6 [> 4 ] 


6.5 [> 4.6 ] 


1.1 [0.7-1.4] 


1.0 [0.7-1.3] 


17: 


X2 


19:43:11.0 


4-23:44:16.1 


10 


3.9 


9.6 [> 6 ] 


11.1 [> 7 ] 


2.1 [1.5-2.8] 


5.9 [4.5-8.0] 




X3 


19:43:09.6 


4-23:43:57.7 


19 


3.5 


10: [> 6 ] 


11.6:[> 7 ] 


2.4 [1.9-3.0] 


2.9 [2.3-3.8] 




X4 


19:43:09.4 


4-23:44:01.0 


17 


2.5 


3.2 [2.2-5.3] 


3.7 [2.6-6.2] 


1.2 [0.8-1.7] 


1.1 [0.9-1.4] 




X5 


19:43:12.5 


4-23:44:21.8 


9 


5.3 


9.3 [> 5 ] 


10.8 [> 5.8 ] 


2.6 [1.6-3.5] 


10: [> 7.2 ] 




X6 


19:43:11.8 


4-23:43:53.0 


12 


3.8 


7.4 [> 5 ] 


8.6 [> 5.8 ] 


1.8 [1.2-2.3] 


3.7 [2.6-5.3] 


18 


X7 


19:43:11.5 


-^23:43:36.5 


7 


3.6 


7.5 [> 4 ] 


8.7 [> 4.6 ] 


0.8 [0.5-1.1] 


2.7 [1.8-3.9] 


12 


X8 


19:43:10.6 


-^23:44:56.8 


12 


2.2 


1.4 [1.1-1.9] 


1.6 [1.3-2.2] 


2.2 [2.0-2.8] 


3.0 [2.5-3.7] 


17 


X9 


19:43:08.6 


4-23:44:24.7 


9 


2.8 


1.2: 


1.4: 




0.3: 


12 


XIO 


19:43:08.9 


4-23:44:36.6 


19 


2.2 


1.5: 


1.7: 




0.8: 


3 


Xll 


19:43:07.7 


4-23:44:49.3 


7 


1.2 


1.0: 


1.2: 




0.2: 




X12 


19:43:14.3 


-^23:45:14. 9 


12 


1.7 


1.3 [1.1-1.7] 


1.5 [1.3-2.0] 




0.6 [0.4-0.8] 


1 


X13 


19:43:13.4 


-^23:45:15.9 


13 


1.2 


1.3 [1.1-1.5] 


1.5 [1.3-1.7] 




0.3 [0.2-0.4] 


2 



3.2.2. Variability 



counts due to statistical fluctuations in a period for which 
two photons are expected, is given by 



We do not see strong intensity flaring in our observa- 
tions (Fig. ^), but sources X2, X3 and X4 show absence 
of detected photons in the time interval between 6 x 10^ 
and 10^ sec, which indicates possible variability in these 
sources. The other sources show a rather constant count 
rate within the 20 ksec of our observation. To quantify 
the possible variabiUty, we performed a statistical test 
of the photon arrival times following the approach out- 
hned by [Preibisch fc Zinnecker (2002| ): given the number 
of detected photons in one source in the total integration 
time, we compute the mean count rate and specify the 
time period in which two photons are expected assuming 
a constant count rate. Then we determine the maximum 
number of photons iV^ax detected in this source within 
any time interval as long as the one derived for the two 
photons. The Poisson probability to find A^max or more 



N, 



k=0 



The probability of variability (pov) is given by pov — 
1 — P (for details see, e.g., Bevington & Robinson 1992). 
The pov values for X3 and X4 are only 68% and 85%, giv- 
ing just marginal evidence for variability within our total 
observing time. Contrasting to that, pov equals 98.3% for 
X2, therefore this source is consistent with X-ray variabil- 
ity. 

3.3. Mass estimates from the near-infrared data 

The 2MASS near-infrared K-band data provide a tool to 
estimate the masses of the infrared sources. Taking the 
observed K-band brightnesses mx at the given K-band 
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extinction (^k ~ 0.11 x A^) we get an estimate of the in- 
trinsic brightness Mk of the source assuming the distance 
d = 2 kpc: 



rriK — Mk = 5 log(c?) 



At 



Different pre-main-sequence tracks are discussed in the 
Hterature (e.g., D'Antona & MazziteUi 1994, Baraffc et al 



199S|, p ^alla fc Stahlcr 1999| , ^iess et al. 200q , p'Antona 
et al. 2pO0D. While some of them differ significantly at the 
low-mass end, they are bet ter comparable going to masses 
> IM© dSiess et al., 20001 ). Most of the calculate d tracks 
do not cover stars more massive than 2.5 M©, but Palla & 



Stahlei (1999) calculate the pre-main-sequence evolution 
to 6Mq and |Siess et al. (2000| ) to TM©. At the high- 
mass end, which we are particularly interested in, both 
calculations agree well. The tracks of intermediate-mass 
stars (> 2M0) are following almost horizontal tracks of 
equal luminosity in the Hertzsprung-Russel diagram, and 
the luminosity is almost independent of the temperature. 
Thus, it is possible to estimate luminosities and by that 
masses from our K-band observations, but the data do not 



nosity of 10'* Lq corresponds in the case of main sequence 
stars to a Bl star (13 Mq). As a cluster of stars is forming 
in IRAS 19410H-2336, the mass distribution of the clus- 
ter members is likely to follow an initial mass function. 
In spite of their larger number, the luminosity-mass rela- 
tion for low-mass stars is lower (L oc m^'* for m < 1 Mq) 
than for high-mass stars (L oc for 1 < m < 30 M0, 
^chatzman fc Praderie 1993| ), and still the major contri- 
bution of the total luminosity stems from the massive clus- 
ter members. For a detailed IMF discussion of the initial 
source sample see ^ridharan et al. (2002 ). As the most 
massive object is likely the central mm source, which is 
not detected in X-ray emission, the X-ray detected sources 
are < 10 Mq. Additionally, the models give the approxi- 
mate luminosity of the sources and we get an estimate of 
the Lx/Lboi ratio (Lx is taken from §3^ and Table |l]). 
The parameters and results are listed in Table 0. 



Table 2. Near-infrared parameters of most likely clus- 
ter sources with 2MASS counterparts. Given are masses (< 
10 M0), bolometric luminosities Lboi, and the ratio Lx/iboi 



give any information about the age of the objects, piess et 
al. (201^) stress that any age determination below 10^ yr 
is highly un certain. In the followin g, we use the tracks 
compiled by Palla fc Stabler (1999), but the resu lts are 
similar when using the tracks of Siess et al. (200C ) . 

Comparing the intrinsic K-band brightnesses A/k with 
the theoretical tracks by Palla fc Stabler (1999), we can de- 
rive mass estimates for the near-infrared objects. We apply 
this procedure for the near-infrared sources corresponding 
to XI, X6, XT and X8, because only those are 2MASS de- 
tections and likely belong to the forming massive clusters 
(j |3.2| ). [- The pro main ocquonco tracka provided by Palla 
& Stal 



src 


niK 




Ak 


Mk 


mass 


ibol 


[10-1 




[mag] 






[mag] 


[Mq] 


[Lq] 


XI 


11.7 


17 


1.7 


-1.5 


> 2 


> 10 


< 70 


X6 


11.0 


18 


1.8 


-2.3 


> 6 






X7 


8.8 


12 


1.2 


-3.9 


> 6 






X8 


12.8 


17 


1.7 


-0.4 


> 3 


> 30 


< 2 



The derived masses are in the intermediate-mass 
regime of Herbig Ae/Be stars, cle arly more massive than 
usual T Tauri or Class I sources. Zinnecker fc Preibisch 



lor (1990 | ) for otaro up to 6 M(.) do not comprioo 

1.6 mag in the K 



sourcc ii intrinoically brighter than 
band. Sourcoo more maooivo than 8 Mm do not have an 



(1994 ) observed X-ray emission from Herbig Ac /Be stars, 
and those sources showed rather soft spectra. [Preibisch 



optically visible pre-main-sequence phase, because they 
start nuclear burn ing already being deep ly embedded in 
their natal cores ( Palla fc Stabler, 1993[ ). Therefore, X6 
and X7 could be more massive than 6 Mq, but it has to 
be noted that the derived model brightnesses Mk neglect 
luminosity due to the disk, i.e., heated dust in the disk, 
and accretion luminosity within the disk and between the 
disk and the protostar. Palla fc Stabler (199E ) stress that 
the understanding of these processes is still inadequate 
for quantitative predictions, but as the near-infrared ex- 
cess in X6, X7 and X8 is not very high (Fig. ||), to our 
j udgm fflnt the additional disk luminosit}" in those sources 



fc Zinnecker (1995 ) speculate that the observations could 
be explained by coronal emission due to a shear dynamo. 
As the spectral signatures of the sources we are study- 
ing here are significantly harder than those found in their 
sample (§ |3.2D , we suggest that the X-ray emission could 
be due to magnetic star-disk interactions as proposed for 
very young Class I sources (Montmerle et al., 200C). In 
that framework, it is possible that the sources observed 
in IRAS 194104-2336 are precursors of the better studied 
Herbig Ae/Be stars. The possible ix/iboi ratio-regime 
(Table ^ for two sources is at the upper boundary con- 
sistent with ratios found in very young low-mass objects 
(~ 2 X 10'"^, e.g., [Feigelson fc Montmerle 1999], [Preibisch fc 



is not ihe dominant effect. This is different for source XI 
which has considerable infrared excess and makes K-band 
bright nesses difficult to be compared with theoretical pro — 
diction s. But as stars below 2 Mq are at least 2 magni 
tudes fainter in the theoretical predictions than the ob- 
served brightness of XI, 2 M0 can still be regarded as a 
tentative lower mass limit for this source. Real upper mass 
limits are given by the total bolometric luminosity of the 
whole cluster of ~ lO''^ Lq (Sridharan et al., 2002). A lumi- 



Zinncckcr 2002 ) , but at the lower end it is also consistent 
with the Lx/-Zvboi ratios found in samples of Herbig Ae/Be 
stars (Zinnecker fc Preibisch, 1994; Gregorio-Hetem et al. J 



1998). It could be argued that the X-ray emission is due 
to low-mass counterparts, but as the derived X-ray lumi- 
nosities are significantly larger than typical values in the 
low-mass regime, most of the X-ray emission stems likely 
from the intermediate-mass sources. 
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3.4. The X-ray sources in detail 

Source XI and the main mm core: The most remarkable 
X-ray source is XI, which does not only show strong 
reddening but also the largest i nfrared excess indicative 



embedded protostellar objects ( Fcigclson fc Montmerle 



1999|)~|As the infrared derived mass is > 2 Mq but the 
X-ray spectrum harder than usual Herbig Ae/Be stars 
( IZinncckcr fc Prcibisch, 1994| ), XI is probably not a nor- 
mal Herbig Ae/Be star but a precursor of such an object. 

Source XI is the X-ray source nearest to the core center 
but not coincident with a mm core. Figure ^ shows that 
the main mm core coincides with a compact cm source 
as well as with H2O and Class 11 CH3OH maser emis- 



sion dMinier et al., 200l| ; [Bcuthcr et al., 2002c) ), while the 
X-ray source is clearly offset and coincides with a near- 
infrared source ~ 6" to the west. Thus, the main mm 
core is likely to contain the youngest and most massive 
object of the cluster, and we do not detect X-ray emis- 
sion from this source, but only from another object nearby 
that might be slightly more evolved. Therefore, we do de- 
tect X-ray emission from pre-main-sequence sources but 
not from the youngest and massive center of the cluster. 
Assuming optically thin dust emission, we estimate the 
core mass of the central source from the 2.6 mm contin- 
uum data to about 80 M0 and the H2 column density to 
w 2 X 10^'' cm~^, corresponding to a visual extinction 
of appro ximately 2000 (for de tails on mm dust calcula- 
tions see Beuther et al. 2002a ). Taking into account this 
visual extinction and assuming different Raymond-Smith 
plasma temperatures for the X-ray emission, we can cal- 
culate a range of upper limits for the X-ray luminosity of 
this source: a 3 keV plasma temperature corresponds to 
an upper limit of Lx < 9 x lO^'' ergs~^, while a temper- 
ature of 10 keV implies Lx < 4 x 10'^'* ergs^^. A lower 
plasma temperature of 1 keV results in an upper limit of 
Lx < 14 X 10'^^ ergs~^. However, it has to be noted that 
the possibility remains that there is no significant X-ray 
emission in the earliest stages of massive star formation. 




11 

E.A. [J2000.0] 



Fig. 4. Zoom into the center of the field presented in Fig. |l| 
Contours and symbols are the same as in Fig. |l| Additionally, 
the black contours on white ground right at the center show the 
cm-peak, the triangle pinpoints the H2O maser position and 



the square the CH3O H maser positions (Minier et al., 2001 
Beuther et al., 2002c|). 



Sources X3 and X4: These sources are both close to the 
same near-infrared source, but accurate astrometry indi- 
cates that X4 rather than X3 is associated with it (Fig. 
^). Additionally, we find an optical counterpart for X4 in 
the Digitized Sky Survey provided by the Space Telescope 
Science Institute. Based on the X-ray spectrum, which is 
not as hard as the spectra of the other cluster members, 
and the optical and infrared data, we cannot distinguish 
whether X4 belongs to the cluster or is an unrelated fore- 
ground object. In contrast, X3 might be associated with 
more diffuse K-band emission (Fig. U). As the hard X-ray 
spectrum of X3 is indicating an early evolutionary stage, 
this source most likely is a very young pre-main-sequence 
object belonging to the cluster under investigation. 



Sources X2 and X5: Source X2 is at the edge of the south- 
ern mm emission and associated with some faint K-band 
emission that might be a reflection nebula from the main 
sources to the south. As the K-band emission is weak, it 
is not detected by 2MASS. Being in the near vicinity of 
a mm core and showing a very hard X-ray spectrum, it is 
possible that X2 might be powered by an embedded pro- 
tostellar object. It has to be noted that this source is the 
only one with observed X-ray variability. Source X5 is also 
too faint in the near-infrared to be detected by 2MASS, 
and it is difficult to decide whether the nearby faint K- 
band emission is a stellar feature or a reflection nebulae. 
But as X5 has the hardest X-ray spectrum of our sam- 
ple, it is very likely a young protostellar object within the 
cluster of IRAS 19410-^2336. 



Sources X6 and X7: These sources also show hard X-ray 
spectra, and both are within the large-scale region of the 
southern massive core as seen in Fig. ^top), but no mm 
core corresponds to any of them. This is reflected in the 
color-color diagram, where they inhabit a reddened region 
without strong infrared excess. The infrared-derived mass 
estimates are not conclusive, but the data - X-ray and 
near-infrared - are consistent with precursors of Herbig 
Ae/Be stars of masses between 5 Mq and 10 Mq. 

Source X8: Source X8 is associated with the strongest mm 
peak of the northern cluster as well as with a near-infrared 
and even an optical counterpart (as found in the Digitized 
Sky Survey). The 2MASS data show that the infrared 
source is reddened and has infrared excess, and we derive 
a lower mass limit of approximately 3.5 Mq. The X-ray 
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spectrum is softer than for the southern cluster members, 
thus it is possible that X8 is in a slightly more evolved 
state of evolution. 



Sources X9 and XIO: The X-ray spectrum of these sources 
is softer, but they are still reddened in the color-color di- 
agram and X9 even shows infrared excess. It is likely that 
they are also young stellar objects, but we are not able to 
determine whether they belong to the cluster or whether 
they are foreground objects. 

Sources Xll, X12 and X13: These sources also have a 
rather soft spectrum, and the color-color diagram locates 
them on or near the unrcddcncd main sequence (Xll is 
too weak to be a 2MASS detection). As they are also spa- 
tially offset from the large-scale mm core, we regard them 
as not associated with the massive star-forming region. 

4. Conclusions 

Hard X-ray emission from a number of point sources is 
detected in a young, massive and embedded star-forming 
cluster in a very early stage of evolution. We did not de- 
tect X-ray emission from the most massive and central 
object (upper limit ix < 9 x 10^"' ergs^^ for kT = 3 kcV) 
but from a few sources in its vicinity. Typical X-ray prop- 
erties of high-mass main-sequence stars {Mi, > 8Mq), 
where the emission originates from internal shocks in the 
radiation-driven stellar winds, are Lx/^boi ~ 10~^, soft 
X-ray spectra with typical temperatures of kT ^0.5 keV 
and very little variability (Bcrghofer et al., 1997). The X- 



ray properties we find for XI to X7 are clearly different, 
they show hard X-ray emission (kT > 3 keV) and most 
of them also near-infrared excess, very similar to those 
observed for extremely young embedded objects (Class I 
protostars), whic h have typical plasma tempe r atures of 
kT 5 - 10 keV ( |Feigelson fc Montmerle, 1999| ; [imanishi 



et al., 2001). Combining infrared data with pre-main- 



sequence evolutionary tracks (|Palla fc Stabler, 1999|) , it 
is possible to estimate the approximate masses of some 
of the hard X-ray sources. Those estimates indicate that 
they are in the intermediate-mass regime of Herbig Ae/Be 
objects. Taking additionally into account the harder X- 
ray spectra compared with other Herbig Ae/Be studies, 
it is hkely that the X-ray sources in IRAS 19410-f 2336 
are even precursors of Herbig Ae/Be stars. The emitted 
X-ray photons with energies mostly above 2 keV indi- 
cate plasma temperatures > 10^ K and X-ray luminosities 
around a few times 10^^ ergs^"'^. The latter values are well 



within the regime of Class I low-mass protostars (Feigelson 
fe Montmerle, 1999), but they are also consistent with 



the results obtained for Herbig Ae/Be stars (Zinnecker fc 
Preibisch, 1994). Thus, some of the objects are probably 



very young intermediate-mass pre-main-sequence sources, 
whereas other sources could also be low-mass Class I or T 
Tauri stars. The emission of one of the sources is consistent 
with X-ray variability. 



In spite of the observation of hard X-ray emission in 
the weak-lined T Tauri star V773 ( |Tsuboi et al., 1998| ), 
where the disk has already been dissipated to a large 
degree, it is unlikely that the hard X-ray spectra ob- 
served in younger class I sources are due to enhanced 
solar-type magnetic activity. Therefore, it is proposed that 
the hard X-ray emission, which is more often observed in 
class I sources than in weak-lined T Tauri stars, is pro- 
duced by magnetic reconnection effects between the pro- 
tostars and their accretion dis ks (pSayashi et al., 19*96 ; 



Feigelson fc Montmerle, 1999| ; [Montmerle et al., 2000| ) 



As the X-ray spectra of the intermediate-mass objects in 
IRAS 19410-1-2336 exhibit very similar signatures to such 
low-mass sources, our results are consistent with disks be- 
ing present in intermediate-mass star formation as well. 

For a better understanding of the nature of the under- 
lying X-ray powering sources much work has to be done 
in the future. Deeper X-ray and near-infrared images will 
help to set stronger constraints on the physical properties 
of the sources: it will be necessary to obtain sensitive X-ray 
spectra to determine better the absorbing A^h column den- 
sities and plasma parameters. It is also of great interest to 
further investigate the properties of the central and deep- 
est embedded object, which means lowering the detection 
limits. Furthermore, the variability of the X-ray sources 
in very young massive star-forming regions is not known 
so far. Therefore, several approaches should be followed 
in the years coming: deep Chandra observations of the 
source of interest will disclose variabilities and faint emis- 
sion of the central object. Additionally, a sample of similar 
sources has to be identified, because only a statistical anal- 
ysis of several young high-mass star-forming regions can 
build a solid picture of the relevant physical processes. As 
high spatial resolution is essential for many of this studies, 
Chandra is a very promising choice. But considering the 
higher sensitivity of XMM-Newton, it might be possible 
to study grating X-ray spectra of the brightest sources of 
the sample of clusters studied then. On the near-infrared 
side, we suggest to get deeper images in the J, H and K 
bands to improve the mass estimates of the X-ray emit- 
ting sources, and near-infrar ed spectroscopy mig ht help 
classifying the types of stars ( Hanson et al., 2002| ). 

To summarize. X-ray studies of young massive star- 
forming regions are just in its infancy, and the next years 
with the space telescopes Chandra and XMM-Newton will 
bring many new insights in that research area. We also like 
to stress that multi-frequency studies over a wide range of 
bands are extremely promising approaches for the under- 
standing of the physical processes forming massive stars. 
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